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ABSTRACT 

Nonlinear redshift-space distortion known as the Fingers-of-God (FoG) effect is a major sys¬ 
tematic uncertainty in redshift-space distortion studies conducted to test gravity models. The 
FoG effect has been usually attributed to the random motion of galaxies inside their clusters. 
When the internal galaxy motion is not well virialized, however, the coherent infalling motion 
toward the cluster center generates the FoG effect. Here we derive an analytical model of the 
satellite velocity distribution due to the infall motion combined with the random motion. We 
show that the velocity distribution becomes far from Maxwellian when the infalling motion 
is dominant. We use simulated subhalo catalogs to find that the contribution of infall motion 
is important to massive subhalos and that the velocity distribution has a top-hat like shape 
as expected from our analytic model. We also study the FoG effect due to infall motion on 
the redshift-space power spectrum. Using simulated mock samples of luminous red galax¬ 
ies constructed from halos and massive subhalos in N-body simulations, we show that the 
redshift-space power spectra can differ from expectations when the infall motion is ignored. 

Key words: cosmology: dark energy - large-scale structure of Universe - galaxies: kinemat¬ 
ics and dynamics 


1 INTRODUCTION 


The peculiar motion of galaxies imprinted on redshift-space clus¬ 
tering provides a good probe of the dynamics of galaxies. The 
coherent motion of galaxies associated with gravitational evolu¬ 
tion squashes the spatial distribution of galaxies along the li ne-of- 
sigh t direction, which is known as the Kaiser effect (IKaiseal 19871 : 
Hamilton 1992j) and provides a unique probe of the growth rate 
jPeacock et al.l200ll ). The growth rate has been measured from the 
results of a variety of galaxy surveys to test g e neral relativity and 


modified gravity (e.g. , lYamtunotoetalj 20081: Guzzoet_aL[ |ZUUS|; 
Ide la Torre et alj|2013l : IBeutler et al. 20141 : Reid et all 2014 ). One 


120081 : 


can expect further precise measurement of the growth rate in a wide 
range of redshifts fro m the results of futu re cosmol ogical surveys 
such as Subaru/PFS 1 Takada et al] I20l4) . Euclid 1 Laureiis et alj 
l201ll) . DESI dLevi et alj|2013l) and WFIRST dSnergel et al.ll20l5tl . 

The internal motions of galaxies in their host halos elongate 
the distribution of galaxies along the line-o f-sig ht direct ion, which 
is called the Fingers-of-God (FoG) effect I Jacksonlll972l) . The FoG 
effect depends on the type of galaxy dZehavi et alj 2005 ). For red 
galaxy samples, such as luminous red galaxies (LRGs), the FoG ef¬ 
fect is influential at scales larger than 10 h ~ 1 Mpc and can introduce 
serio us systematics in the measur e ment of the cosmic growth rate 
(e.g., iHikage & Yamamotol 120131 : IBeutler et al.l l2014l : lReid et alj 
120141) . Clustering anisotropy can generally be expanded with a 
series of multipole power spectra Pi(k) (Yamamoto et ah 2006). 
The Kaiser effect mainly generates quadrupole anisotropy ' Kaiser] 


Il987l) . Meanwhile, the FoG effect generates higher multipole 
anisotropy, such as hexadecapole (l = 4) and tetra-hexadecapole 
(l = 6) components, whe re the Kaiser effect is subdominant 
dHikage & Yamamotd [2013). These multipoles are useful in elim¬ 
inating the FoG uncertainty and provide kinematic in formation of 
satellite galaxies faikage)l2014l : lKanemaru et al]|2015t) . 

The FoG effect is related to the kinematics of central and 
satellite galaxies inside their host halos. The satellite dynamics 
can be altered by diff erent phys ical processes such as dynami- 
cal friction (Chandrasekhar 1943 }), tidal stripping/ disruption (e.g., 
iBovlan-Kolchin et aPl2008l : Wetzel & WhitelioTol) . satellite merg- 
ing, and hydrodynamical drag, such as ram pressure l lGunn & Gotti 
1 1972n . Such complicated processes of satellite kinematics intro¬ 
duce systematic uncertainties in the determination of the dynam¬ 
ical mass of g alaxy clusters through satellite velocity dispersion 
(e.g.. IWu et al.ll2013l) . The velocity distribution of satellites has 
been studied by conducting various numerical simulations (e.g. , 
Ghignaet alj |200d : iDiemand et al Jl2004l : iFaltenbacher et al]|2005l : 

Wuetalluom A non-zero velocity of central galaxies rela- 

tive to the host halo has been found in several studies (e.g., 
van den Bosch et al . 200 llGuo etaPl2015l) . 

The FoG effect has usually been attributed to the random mo¬ 
tion of galaxies inside clusters or halos. Meanwhile, coherent infall 
motion onto the h alo mass ce nter also generates the FoG effect, as 
illustrated by iHamiltonl d 1992h . The infall region of clusters form¬ 
ing a trumpet-shaped pattern has been observed for a number of 
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Figure 1. Illustration of the line-of-sight velocity distribution dlos due to random motion (left), coherent infall motion (right) and their combination (center). 
The velocity distribution due to coherent infall motion has a top-hat shape, while the random motion produces a Gaussian velocity distribution. 


clusters and provides an important probe of the clu ster mass pro¬ 
file insensitive to the details of galaxy formation i Regos_&_Gellei 
19891 : iDiaferio & Gel let! 1 1 9971 : iRines et al]|2003l : Zu & Weinbera 
2013i) . In this letter, we study the FoG effect due to the infall mo¬ 
tion of satellites on the redshift-space power spectrum. We derive 
a theoretical model of the line-of-sight velocity distribution due to 
the random and infall motions of satellites inside their halos. We 
use simulated subhalo catalogs to test our model expectations and 
to study the effect on the redshift-space power spectra. Flere we 
focus on LRGs that have been widely used in cosmological stud- 
ies (e.g jReid et all2()Iol : lBeutler et al.l2014l : lReid et al.l2014l) . Al¬ 
though the satellite fraction of LRGs is only about 6%, the effects 
of the satellite FoG e ffect on the growth rate me asurement have 
been found to be large iHikage & Yamamotol20131) . Understanding 
the behavior of satellite motion is important in precise cosmologi¬ 
cal studies using redshift-space power spectra of LRG samples. 

This letter is organized as follows. In section [2] we model 
satellite internal motion by decomposing the infall and randomly 
rotating components. The details of our simulations are explained 
in section[3] In section[4] we present results, comparing the satellite 
velocity distribution with the theoretical modeling derived in sec¬ 
tion [2] and evaluate the FoG effects of the infall velocity and ran¬ 
dom velocity on redshift-space galaxy clustering. Section [5] sum¬ 
marizes the letter and presents conclusions. 

2 SATELLITE MOTION INSIDE HOST HALOS 

The internal satellite velocity with respect to the host halo bulk 
velocity v sa t = V aa t —Vh a i 0 can be decomposed into a component 
infalling onto the halo center and a tangential component: 

Vsat — ( (tti n f) F 6inf)Gr- "F Ct an ,0G$ F £tan, 4 >^</}i (1) 

where the mean infall velocity (ui n f) is generally non-zero and de¬ 
fined to have a positive sign in the direction toward the halo center. 
The average velocity dispersions in the infall and tangential direc¬ 
tions are defined as o^ inf = (e? nf ) and cr^ tan = (e? an ,s) = 
( e t an ,^) respectively, and depend on the host halo mass M. The 
line-of-sight (one-dimensional) component of the internal satellite 
velocity ulos = v sat • e L os becomes 

vlos = (- (ttinf) F einf)/x F et an ,s(l - ^ 2 ) 1/2 , (2) 


the infall and tangential velocities are given by independent Gaus¬ 
sian distributions, the line-of-sight velocity distribution of satellites 
in the direction of /j becomes a Gaussian distribution with a mean 
of —/r (ttinf) and variance of 

Ov,m(M; AT) = ft 2 OV, inf(AT) F (1 - F 2 )^,tan (AT). (3) 


The line-of-sight velocity distribution averaged over ^ is given by 
d/r- 


/v(vlos;AT) = - 


x exp 


(27r) 1 / 2 CT v , M (/r; AT) 

('bLOS F H (finf) (AT)) 2 
2o-v, M (ft; AT) 


(4) 


The velocity dispersion of satellites is the second moment of the 
equation 0 and becomes 

°v,LOs(AT) = -((ttinf) 2 F Ov,inf + 2cr^ tan ). (5) 

When the internal motion of galaxies is well virialized, the 
mean infall velocity is negligible compared with the random mo¬ 
tion and the internal velocity distribution is isotropic, i.e., cr v ,i n f = 
ov,tan = ovir, where ovir is the virial velocity dispersion. Equation 
0 simply becomes the Gaussian distribution 


fv(v los; AT) = 


(2tt) 1 /2 c 


- exp 


2cr v,vir 


( 6 ) 


which is illustrated in the left panel of Figure 0 Meanwhile, when 
the mean infall velocity is dominant ((fmf) (J v ,i n f, rr V! t an ), 
equation 0 becomes a top-hat distribution: 


/u(vlos; AT) 


0.5(t?inf) (IttLOsI ^5 (t^inf)) 

0 (I'ULOsj > (tFif)) 


(7) 


where the line-of-sight velocity dispersion becomes ct v ,los = 
(finf)/\/3 (see the right panel in Figure 0 . In reality, the infall 
motion and random motion combine and the line-of-sight velocity 
distribution is then described by the sum of Gaussian distributions 
with different mean values of fi (rii n f) in the range of ^ from —1 to 
1. The shape of the distribution then becomes a smoothed top-hat 
function like that seen in the center panel of Figure 0 


3 SIMULATIONS 


where /r is the cosine of the angle between the line of sight and the To validate our modeling of the line-of-sight velocity distribution 

direction of infall toward the halo center; i.e., /r = e r • eLos- When derived in the previous section, we construct subhalo catalogs from 
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Figure 2. Halo mass dependence of the line-of-sight (one-dimensional) velocity dispersion <x v lqs (black circles), radial and tangential velocity dispersions 
tr v ,inf (yellow square) and <x v ,tan (red triangle) and mean infall velocity (ui n f) (blue crosses) relative to the virial velocity dispersion rr v ir- 


N-body simulations. We make 10 realizations using GADGET-2 
public code i jspringefcoosi) . starting from an initial redshift of 80 
and a 2LPT initial condition ICrocce et al .1120061) . The side length 
of the simulation box is 600 ft , -1 Mpc and the number of particles is 
800 3 , which corresponds to a particle mass of 2.8 x 10 lo /i _1 M©. 
The softening length is set to 10/U 1 kpc. We_compute the initial 
power spectrum using CAMB software (Lewis_et al.l l2000l) by as¬ 
suming a flat ACDM cosmology with the parameters fl m = 0.273, 
fib = 0.046, h = 0.704, n s = 0.963, r = 0.089, and os = 0.809. 
We use snapshot data at z = 0.3 roughly corresponding to the 
mean redshift of the LRG s ample of the Sloan Digital Sky Sur¬ 
vey (SDSS) feisenstein et alj|200ll) . Halos are identified using the 
friends-of-friends algorithm with a linking length b = 0.2 and 
setting the minimum number of particles fo r halos to 32. Subha¬ 
los are identified using the SUBFIND code ISnringel et al.ll200ll) 
with the minimum number of particles set to 20. We take the z-axis 
as the line-of-sight direction to add the peculiar velocity effect as 
r —> r + (1 + z)v z /aH(z). 

We describe the relationship between h alos and galaxi es us- 
ing the halo occupation distribution (HOD) l lKravtsov et al]|2004c 
IZehavi et all2005l : IWhite e t all 201 ll). We use a conventional HOD 
form with five parameters (IZheng et alJi2005h : 


(Alcan (M)) 


(Asat(M)) 


1 1 + erf / loglo(M) ^ log 10 (Minin) 

2 V CTlog M 



M - Af cut 

Mi 



(?) 

(9) 


where erf (a:) is th e error function. We fix the HOD values for an 
SDSS LRG as jReid & Spergell2009l) = 5.7x lO 13 fc _1 M 0 , 
a iosM = 0.7, M cut = 3.5 x lO^ft^M©, Mi = 3.5 x 
10 14 /i _1 Mq, and a = 1. The position of the central galaxy is 
assigned to be the potential minimum of the halo and the veloc¬ 
ity is assigned as the mean of all dark matter particles inside cen¬ 
tral subhalos. We use three different types of tracers for satellite 
galaxies; one type is randomly selected dark matter particles and 
the others are subhalos with mass relative to the host halo mass 
/ su b = Msub/Mh a io of <0.1 and ^ 0.1 (hereafter referred to as 
“light subhalos” and “heavy subhalos” respectively). 

4 RESULTS 

4.1 Satellite velocity distribution 


We use simulated subhalo catalogs and decompose the internal ve¬ 
locity of satellites relative to the host halo velocity into infall and 


tangential velocity components. Figure [2] shows the results of the 
host halo mass dependence of each velocity component: the line- 
of-sight velocity dispersion <j v ,los, mean infall velocity (ui n f) and 
the infall and tangential velocity dispersion cr v ,mf and ov,tan- The 
panels show the result when satellites are represented using dark 
matter particles, light subhalos, and heavy subhalos from left to 
right. Each velocity component is normalized with the virial veloc¬ 
ity dispersion cr v i r for the corresponding host halo mass. 

When satellites are represented by dark matter particles, the 
mean infall velocity is low and the line-of-sight velocity dispersion 
is nearly equal to the virial velocity dispersion. This indicates that 
the motion of dark matter particles is well virialized. Meanwhile, 
the motion of subhalos is quite different from that of dark matter 
particles. The mean infall velocity of subhalos is much higher than 
that of the dark matter particles and it is comparable to the virial 
velocity for massive subhalos. The tangential velocity reduces as 
the mean infall velocity increases. This indicates that the dynami¬ 
cal friction becomes efficient for heavier subhalos and the orbital 
motion slows down. The velocity bias of the line-of-sight velocity 
dispersion varies depending on the host halo mass and the subhalo 
mass. Our result shows that the line-of-sight velocity dispersion is 
smaller than the virial velocity dispersion. The result may change 
quantitatively when the subhalo finders are different. The inner fast 
rotating subhalos are difficult to identify using the SUBFIND al¬ 
gorithm and the average velocity of satel lite subhalos is thus lower 
than that in previous works (e.g jDiemand et ai]|2004l) . The quali¬ 
tative result that the infall motion is dominant in heavier subhalos 
is consistent with the results obtained previously. 

Figure [3] shows the line-of-sight velocity distribution of satel¬ 
lites represented by dark matter particles (left), light subhalos (cen¬ 
ter) and heavy subhalos (right) for different bins of host halo mass. 
The velocity distribution for dark matter satellites agrees with a 
Gaussian distribution with virial velocity dispersion for each bin 
of the halo mass. For satellite subhalos, the mean infall velocity 
is comparable to or larger than the tangential velocity. The result¬ 
ing shape of the line-of-sight velocity distribution does not become 
a simple Gaussian but has a more top-hat-like shape, particularly 
for heavy su bhalos. This feature is consistent with previous results 
obtained by Diemand et ali J20041) . who pointed out that the sub¬ 
halo velocity distribution becomes non-Maxwellian. We compare 
the simulation results with our modeling using the simulated values 
of the mean infall velocity and velocity dispersions. For heavy sub¬ 
halos (the right panel of Figure[3]l, our modeling of the line-of-sight 
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Figure 3. Distribution function of the line-of-sight velocity -ui os distribution of satellites for three different ranges of halo mass from 10 14 /i —1 M© to 
8 x 10 14 /i _ 1 Mq. Satellites are represented with dark matter particles (left), light subhalos (center) and heavy subhalos (right). Reference lines are a Gaussian 
distribution of the one-dimensional velocity dispersion (eq. [6), dashed) and our velocity model including both the infall and random motions (eq. (4), solid). 
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Figure 4. Comparison of the FoG damping function T>(k y) (eq. HOl ) 
when the infall motion or random motion is dominant. For comparison, 
we plot the Lorentzian form of FoG damping. The velocity dispersion 
<r v LOS/(°0 i s set to 1 Mpc. 


Figure 5. Mock LRG power spectra Pi simulated with the halo model de¬ 
scription. Satellites are represented by massive subhalos (where the average 
ratio of the subhalo mass to the host halo mass is 0.06) and the satellite 
fraction is about 6.4%. For comparison, we plot the power spectra where 
the velocity of satellite subhalos is rotated randomly. 


velocity distribution given by equation 0 describes the simulation 
results better than a simple Gaussian distribution (eq. t6l). 


4.2 FoG effect on the redshift-space power spectrum 


Satellite motion inside halos generates non-linear redshift-space 
distortion, often referred as the FoG effect. The FoG effect on 
multipole components of the redshift-space power spectr um is 
formul ated em ploying the halo model (e.g., Hikageet II] 1 20 if 
IHikage & Ya mamo tol l2013l ). FoG damping due to internal satellite 
motion of the halo with mass M is given by Fourier transforming 
the line-of-sight velocity distribution of satellites: 


V(kf,M)= / dvfv(v,M)exp(—ik\\v) 


~ exp 


kfOv.sat - («inf) 2 /3) 


sin(fe|| (ttinf) ) 
k\\ (U in f) 


( 10 ) 


where k = (1 + z)k y /H(z) and fc|| is the line-of-sight component 
of the wavevector k. In the second line, we make the approximation 


that ov^inf ~ (Jv,tan- Figure[4]compares the FoG damping function 
in different cases. When the mean infall velocity is equal to zero, 
FoG damping has a simple Gaussian form. When the infall veloc¬ 
ity is dominant, the velocity distribution has a top-hat form and its 
Fourier-transform is a sine function, as shown by the dashed line in 
Figure |4] FoG damping due to the infall motion has a shape quite 
different from the commonly used Gaussian or Lorentzian form and 
affects the small-scale features of redshift-space clustering. 

Figure [5] shows the multipole power spectra k 2 Pi(k) (l = 
0, 2, 4 from left to right) of simulated samples assuming the FIOD 
of SDSS LRGs. Satellites in the samples are represented by heavy 
subhalos, which is a good assumption because LRGs are mas¬ 
sive galaxies and thus should have massive host subhalos. The av¬ 
erage ratio of the subhalo mass to the host halo mass becomes 
6.4%. In this sample, the mean infall velocity is comparable to 
the virial velocity of satellite-hosting halos and the values of each 
velocity component of satellite subhalos are (unit) = 0.79cr v i r ; 
Cv,mf = 0.95cr v i r ; (T v ,tan = 0.86a v i r - To see the effect of co¬ 
herent infall motion, we also plot the power spectra of the same 
mock samples where the direction of the internal satellite velocity 
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v sa t is randomly rotated to erase the coherent infall motion while 
the absolute value |v sa t| is the same. The error bars represent the 
1-sigma scatter of the power spectra Pi for 10 LRG mock sam¬ 
ples divided by the square root of 10. Assuming that the sample 
variance is proportional to the survey volume, the plotted error cor¬ 
responds to the statistical error for the survey with the sample vol¬ 
ume of 2.16fr -3 Gpc 3 (= 10 x (600/i _1 Mpc) 3 ), wh ich is roughly 
the sa me volume as for the BOSS/CMASS sample dBeutler et alj 
l20l4) . The effect of the infalling motion of satellites becomes ap¬ 
preciable on smaller scales (k it 0.2/iMpc -1 ) because the internal 
velocity distribution is sensitive to the small-scale power spectrum 
via the FoG effect. We evaluate the effect of the infall FoG on the 
overall shape of multipole power spectra as chi-square values: 




2 


E E 


l i 


- Pf /o infall (fc)] 2 


a 


2 

Pi 


(ki) 


(ID 


where op, is the lcr error of Pi for the survey with 2.16 h~ 3 Gpc 3 
volume, and fc max is the maximum value of k for the calcu¬ 
lation. The binning width of k for the calculation is set to be 
0.04/iMpc _1 . The chi-square value increases at larger fc max : \ 2 = 
10 at fc max = 0.21/iMpc _1 , X' 2 = 89 at &max = 0.29/iMpc -1 , 
and x 2 = 190 at fc max < 0.41/iMpc” 1 . Our result indicates 
that the proper treatment of the infall FoG would be necessary for 
the precise modeling of redshift-space power spectra for massive 
galaxies such as LRGs at k it 0.2/iMpc _1 . 


5 SUMMARY AND CONCLUSIONS 

We investigated satellite motion inside host halos and the resulting 
FoG effect on the power spectrum. When the satellite motion is not 
well virialized, the internal satellite motion is not random but has 
coherent infalling flow onto the halo center. We used subhalo cata¬ 
logs to find that the infall motion becomes important when the sub¬ 
halo mass relative to the host halo mass increases and the line-of- 
sight velocity distribution then deviates from Maxwellian and has 
a flat top-hat form. We derived a theoretical model of the satellite 
velocity including both infall and random motions. We found that 
our model well describes the simulated subhalo velocity distribu¬ 
tion. We also investigated how the difference in the velocity struc¬ 
ture between random motion and infall motion affects the shape of 
redshift-space power spectra. We showed that the effect of infall 
motion on the overall shape of an LRG-like power spectrum is ap¬ 
preciable on the scale k > 0.2 h per Mpc for the current survey 
volume, such as the BOSS/CMASS sample volume. 

The effect of infall motion may be important in precision 
cosmological studies using redshift-space galaxy clustering in fu¬ 
ture surveys such as Subaru/PFS, DESI, Euclid, and WFIRST. We 
found that the contribution of infall motion depends on the mass of 
subhalos hosting satellites. Galaxy-galaxy lensing provides the in¬ 
formation of subhalo mass hosting satellite galaxies (e.g jLi et alj 
l2014ll . It would be interesting to see how the galaxy-galaxy lensing 
measurement from the upcoming imaging and redshift surveys can 
remove the uncertainty of the infall motion effect. We leave this as 
a topic of research to be addressed in the near future. 
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